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ABSTRACT 

We search for the dark flow induced diffuse kinetic Sunyaev Zel'dovich (kSZ) effect through CMB- 
galaxy cross correlation. Such angular correlation is anisotropic, with a unique cos(6*df) angular 
dependence and hence can be distinguished from other components. Here, #df is the angle between 
the opposite dark flow direction and the direction of the sky where the correlation is measured. We 
analyze the KIAS-VAGC galaxy catalog of SDSS-DR7 and the WMAP seven-year temperature maps, 
applying an unbiased optimal weighting scheme to eliminate any statistically isotropic components 
and to enhance the dark flow detection signal. Non-zero weighted cross correlations are detected 
at 3.5a for the redshift bin z < 0.1 and at 3cr for the bin 0.1 < z < 0.2, implying the existence of 
statistically anisotropic components in CMB. However, further analysis does not support the dark flow 
explanation. The observed directional dependence deviates from the oc cos #df relation expected, and 
hence can not be explained by the presence of a single dark flow, and if the observed cross correlation 
is generated by the dark flow induced kSZ effect, the velocity would be too high (> 6000 km/s). 
We report this work as the first attempt to search for dark flow through weighted CMB-galaxy cross 
correlation and to draw the attention on the sources of the detected anomalous CMB-galaxy cross 
correlation. 

Subject headings: Dark flow - Kinetic Sunyaev- Zel'dovich - cross correlation - galaxy and CMB 



1. INTRODUCTION 

The peculiar velocity field of galaxies provides an im- 
portant and robust way to understand the matter distri- 
bution of the universe. In the standard cosmology, the 
peculiar velocities arise from matter density perturba- 
tion and are unbiased tracer of the background matter 
distribution. It has been used to constrain the amplitude 
of matter power spectrum or the matter density frac- 
tion f2 m , and it could be used to detect possible sources 
of gravitational field which can not be observed directly 
with existing galaxy surveys. 

The peculiar velocity field of the nearby universe has 
been explored extensively using variant velocity trac- 
ers. On small scales (< 20 ft,~ 1 Mpc) the peculiar 
velocities of ga laxies are consistent with the standard 
ACDM model dPike k Hudson! [20051: lAbate k Erdogdul 
120091 [Juszkiewicz et al.ll201Q[ ). However, on larger scales, 
various groups have reported bulk flows with unexpected 
amplitude. (I) On intermediate scales extending to 
about 100 ft, _1 Mpc, there is mounting evidence of a 
bulk flow with amplitude ~ 400 k ms" 1 using a compo- 
sition of current velocity catalogs (iWatkins et al.l [20091 : 
iFeldman et al.ll201Ct IMacaulav et alj|2011ft . which is un- 
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likely within the frame of standard ACDM model, despite 
the fact that there are some authors who claimed nor- 
mal bulk ve locities on these scales using seperate veloc- 
ity catalog (iNusser k Davisl 120111: iTurnbull et~aLll2012t 
iMa k Scottll2012D . (2) On Gpc scales extending toz~ 



0.3, an even larg er bulk flow of about 1000 km s wa s re- 
ported recently (jKashlinskv et al.l 120081. 120091 120101) . by 
measuring the kinetic Sunyaev-Zel'dovich (kSZ) effect of 
galaxy clusters in WMAP. Whether or no t this dark flow 
exists is stil l under intensive detat e (c.f. lOsborne et al.l 
(12011ft and iModv k Haiianl (|2012ft for negative result, 
and lLi et al.l (|2012l ) for peculiar velocities in ACDM sim- 
ulation). If confirmed, this large bulk flow, dubbed "dark 
flow", would severely challenge the standard model of 
cosmology and structure formation. 

It is therefore valuable to perform independent tests 
of the dark flow scenario beyond the cluster kSZ effect. 
In this paper we make an independent observational 
check of the dark flow scenario by cross-correlating 
the SDSS-DR7 galaxies with the WMAP 7 year CMB 
data. The galaxy-C MB cross correlati o n has been 
measu red extensively dFosalba et all 120031: iCabre et aT 
2006aibt rGiannantonio et all 12000. I2008t IHo et al 
2008; iHirata et al.l 120081: | L6pez-Corredoira et al 
20ld iRaccanelli et all [20081: iSawangwit et al 



20101 iFrancis k PeacocM 120101: iBielbv et a' . 
Hernandez-Monteagudol 12010ft . mainly to detect the 



integrated Sachs-Wolfe (ISW) effect. Similar analysis 
can be carried out to probe da rk flow, under a specially 
designed weighing scheme. As Zhang (2010) (hereafter 
Z10) proposed, a large dark flow will induce significant 
small scale anisotropics on the CMB temperature 
through the kSZ effect, due to fluctuations in the free 
electron number density. On each small patch of the 
sky, the kSZ effect induced by the clustering pattern of 



2 



this dark flow traces the large scale structure closely. 
However, the overall amplitude is modulated and has a 
dipole pattern across the sky, with the positive maxima 
located at the direction opposite to the dark flow, and 
the negative minima located at the direction of the 
dark flow. This statistically anisotropic distribution is 
distinctively different from all cosmic sources such as 
the primary CMB, the ISW effect, the thermal Sunyaev 
Zel'dovich (tSZ) effect, and the cosmic radio and infrared 
backgrounds. For this reason, one can apply a weighting 
function of appropriate directional dependence to elimi- 
nate the (statistically) isotropic components and isolate 
the kSZ component in CMB-galaxy cross correlations. 
This is a major difference between our work and other 
existing works on galaxy-CMB cross correlation. The 
proposed method can be extend to any redshift to 
measure coherent large scale flow against CMB, if any. 
It is comp l ement ary to the recent measurement by 
lHand et all (|2012D . which measures the mean pairwise 
momentum of clusters through the kSZ effect of clusters 
and hence probe the peculiar motion at smaller scales. 

In §2 we review the theoretical prediction of the dark 
flow induced kSZ effect. The CMB temperature maps 
and galaxy catalog we used are described in §3, and the 
data analysis is done in §4. Results are shown in §5, 
along with investigations on possible systematics. We 
conclude the paper in §6, with a brief summary. 

2. THE KSZ-GALAXY CROSS CORRELATION 

Dark flow, if exists, will induce CMB temperature fluc- 
tuations at all angular scales through the modulation of 
the inhomogeneous electron distribution on the uniform 
dark flow (Z10). The induced kSZ effect is given bj0 
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where Tcmb = 2.725K is the CMB temperature at 
present, Vdf is the dark flow velocity vector, with Vdf 
as its amplitude and #df the angle between the line of 
sight f and the opposite dark flow direction. The effective 
projection function for the kSZ signal is defined as 



zf^L , (2) 



where Xe is the ionization fraction, <tt is the Thomson 
cross section, and Qb is baryon density parameter. The 
integration in Eq. ((T|) is over the dark flow region along 
the line of sight. The first term in the last expression is a 
pure dipole term, which has no correlation with the large 
scale structure. For this reason, we will neglect this term 
hereafter. The second term represents the kSZ effect 
arising from modulation of the inhomogeneous electron 
distribution (S e ) on the dark flow through the inverse 

8 The original paper (Z10) missed a minus sign. However, nu- 
merical results presented in Z10 is unaffected. 



Compton scattering. Clearly, this component traces the 
large scale structure and is tightly correlated with galaxy 
distribution. Hence it is promising to search for dark 
flow by cross correlating CMB temperature maps with 
galaxy distribution over the redshift range of dark flow. 
The induced cross correlation is 



w Tg {9) = (AT ksz (n)^(n ))| 



(3) 



Here, 5g is the galaxy surface number density fluctua- 
tion, 9 is the angular separation between the two lines 
of sight, with cos 9 = h ■ n . w T9 can be measured over 
either the whole survey area, or in each sub-region of 
it. We denoted this area average by a subscript "S". 6 
should not be confused with 6*df, which is the angular 
separation between the opposite dark flow direction and 
the line of sight. 

The kSZ-galaxy cross correlations induced by dark 
flow have a number of signatures by which it can be 
distinguished from the cross correlations arising from 
other CMB components, such as the thermal Sunyaev 
Zel'dovich effect, the integrated Sachs- Wolfe effect, and 
the foreground contaminations. Below we discuss these 
one by one: 

• w T9 is statistically anisotropic, with a character- 
istic cos(#df) angular dependence. On the con- 
trary, cross correlations induced by the thermal 
SZ effect, the ISW effect and extragalactic fore- 
grounds are statistically isotropic. One can search 
for such signature by splitting the survey sky into 
sub-regions according to the values of cos(6*df) and 
measure w Tg in each sub-regions. Alternatively, 
one can choose an appropriate weighting function 
W{n) such that the weighted cross correlation 



(9) ee (AT^(h)W(h)S^(n )}| 
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vanishes for any isotropic CMB components. This 
requires (W)s = 0. Here we denote this weighted 
cross correlation with a subscript "W" . The opti- 
mal weighting function minimizing the cross corre- 
lation measurement error is3 



W(n) = cos#df — (cos$df)s 



(5) 



The weighing scheme eliminates statistically 
isotropic components. Under this weighing scheme, 
the dark flow induced kSZ survives and the sig- 
nal is optimized. However, anisotropic compo- 
nents of more complicated angular dependence 
(e.g. cos(7j0df) {n = 2, 3, • • •)) also survive under 
the weighting. So a non-vanishing w^y means the 
existence of anisotropic components, but not neces- 
sarily means the existence of dark flow. To further 
test it, we will compare the measured w Tg as a 
function of 6*df and check if it follows the cos(6>df) 
dependence. 

• w Tg should be frequency independent, since kSZ is 
frequency independent. 

9 Z10 proves that the optimal weighting function satisfies the 
equation W = (W 2 )$ [cos 6*df — (cos #rjp),g]/{cos Sop W)s- Here 
we show that the solution to the above equation is W = cos £>df — 
(cos9df)Si u P to an arbitrary prefactor of no 6>df dependence. 
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w w(@) ( anc ^ w T9 (0)) should have similar 8 depen- 
dence as that of the galaxy angular correlation 
function w" . The two cross correlations are re- 
lated by the following order of magnitude approxi- 
mation, 
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Here, fw — (cos(0df)W' (n)) s , r e is the Thompson 
optical depth over the dark flow redshift range, b e 
is the bias of electron number distribution with re- 
spect to matter distribution and b g is the galaxy 
bias. For the dark flo w direction (Zdf = 282° , 
&df = 22°) reported by iKashlinsky et all ((2010), 
fw ~ 0.09 over the overlapping area of KIAS- 
VAGC and WMAP7. A complexity is that b e 
could have non-negligible scale dependence over 
relevant scales due to various gastrophysics such 
as AGN and supernova feedback. So the above 
relation only serve as a guideline. In the follow- 
ing, we use a fitting formula of gas window fr om 
hydrodynamical simulation of Uing et al.l (|2006T ) to 
account for b e . We define the scaling factor as 
ratio of observed galaxy corr elation function t o 
the halofit model prediction ([Smith et al.1 [2003). 
R{9) = w B l s /Wf^ alo f it , and the galaxy bias is sim- 
ply, b g = vR- The resulted kSZ signal induced by 
a single dark flow are shown in Figure [3] as dashed 
lines. 

Hereafter we will search for the dark flow through such 
signatures in the WMAP-SDSS cross correlation. 

3. DATA 

Our cross correlation measurements are based on the 
KIAS-V AGC galaxy catalo g of SDSS DR7 main galaxy 
survey (|Choi et alj 12010a!) and the CMB temperature 
maps from the WMAP seven-year data release. 

3.1. KIAS-VAGC galaxy catalog 

The main source of KIAS-VAGC is the New York Uni- 
versity Value-Added galaxy catalog (NY U-VAGC) large 
scale structure (LSS) sample "brvoidO" dBlanton et al.l 
2005), which has r-band magnitude limit of 17.6 rather 
than SDSS limit 17.7 to maintain a homogeneous sur- 
vey depth and a simp le survey boundary. All magni- 
tudes are k-corrected ([Blanton et al.ll2003al ) and evolved 
to rest frame magnitudes at redshift z=0.1 using an 
updated versi o n of t he evolving luminosity model of 
iBlanton et al.l (|2003bl ). The angular selection function 
is characterized carefully for each sector on the sky. The 
fiber collision results in about 7% of targeted galaxies 
not having measured spectra. These galaxies are as- 
signed the redshifts of their nearest neighbors with avail- 
able spectra. This method has been shown to be rea- 
sonably accurate fo r large scale structure measurement 
(jZehavi et all 12005). Additionally, about 10k redshifts 
are borrowed from other spectroscopy gala xy surveys, 
such as the updated Zwicky Catalog (UZC) (Fal co et all 
11999ft . the I RAS Point Source C atalog Redshift Sur- 
vey (PCSZ) (jSaunders et al.|[2000h . the Third Reference 
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KIAS-VAGC NGC (cross WMAP7) 
Completeness: 0.67 
covered area: 6865 degree 
Effective area: 6518 degree 2 
Galaxy number: 514207 
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Fig. 1. — The galaxy number density with redshift in our home 
catalog. 




0.3 

0(degree) 



0.5 0.6 0.7 



Fig. 2. — Auto-correlation functions of galaxies for z < 0.1 (solid 
line) and 0.1 < z < 0.2 (dashed line). 



Cata log of Bright Galaxies (RC3) (jde Vaucouleurs et all 
119951) . and t he Two Degree field Galaxy Redshift Survey 
(2dF-GRS) (IColless et all 12001 . The resulted catalog 
contains 593,514 redshifts and covers 2.56 sr of the sky. 
In high density regions, KIAS-VAGC has high complete- 
ness. More details about the KIAS-VAGC catalo g are 
given in IChoi et all (12010bft and lChoi et all (12010aft . 

In this work we use the major part of KIAS-VAGC, 
the north Galactic cap (NGC), which covers an effective 
area of 6518 square degree. The sky area overlapping 
with the WMAP7 temperature maps out of the KQ75 
mask contains 514207 galaxies. The number density of 
the resulted galaxy catalog is shown in Figure [TJ To 
study the radial depth over which the dark flow extends, 
we split the galaxies into two redshift bins: z < 0.1 and 
0.1 < z < 0.2. 

3.2. CMB Temperature Maps 

The WMAP seven-year dat a products are pu blicly 
available in the Healpix format (jGorski et alj |2005) . The 
WMAP observed in 5 frequency bands, i.e. the K band 
at 23 GHz, Ka band at 31 GHz, Q band at 41 GHz, V 
band at 61 GHz, and W band at 94 GHz. The first two 
low frequency bands are used to extract Galactic fore- 
grounds. The 3 high frequency bands Q, V and W (with 



4 



8 channels totally) has subdominant foreground contam- 
inations and relatively higher resolutions, suitable for 
testing frequency-dependent contaminations to the small 
scale kSZ-galaxy cross correlations. In fact we do not 
see any considerable difference in the cross correlation 
results and the associated measurement errors of these 
three bands. We do not use the foreground-free ILC ma p 
provided by the WMAP team (jDelabrouille et al.ll2009h . 
due to its low angular resolution (with an FWHM about 
1°). We use the CMB temperature maps with resolution 

9 (res = 9) in Hcalpix format, corresponding to an angu- 
lar resolution of 7 arcmins. This results in 3145728 pixels 
with 47.21 square arcmin each. We use a combination of 
the "Extended temperature analysi s mask" (KQ75) and 
the "Point Source Catalog Mask" (jGold et all 120111) to 
exclude unobserved or foreground contaminated regions. 
About 70% sky regions are available after these mask- 
ings. 

To explore the systematics coming from intrinsic CMB 
fluctuations, and to test the robustness of the Jackknifc 
error estimation method in this analysis, we generated 
1000 CMB temperature simulations with the "Synfast" 
code provided by WMAP. These simulated maps are 
drawn with the best-fit CMB power spectrum, then we 
add thermal noise to each map. The thermal noise are 
Gaussian distributed with a variance inversely propor- 
tional to the observation times of each pixel. The ther- 
mal noise in different pixels are uncorrelated. 

Due to large errors in the cross correlation measure- 
ments, we do not use the cross-correlated mock samples 
for galaxy and CMB. It has been shown that the esti- 
mated errors in the cross correlation measurements us- 
ing the CMB simulation maps are consistent with those 
using cross correla ted mock samples with high precision 
(jCabre et al.ll2007t ). 

4. CROSS CORRELATION AND ERROR ESTIMATION 
The galaxy overdensity in the i-th pixel is = 
n f,i/™g,i~ 1' where is the galaxy surface number den- 
sity and fig f is the expected mean galaxy number density. 
The CMB temperature anisotropy is ATj —Ti — Tqmb- 
The cross correlation function at angular separation 9 is 
given by 



w T9 (0) = ^ 



y.,.j fg,ih,j 



(7) 



Here f g ^ is the galaxy selection function or completeness 
and fx,i is the fraction of pixel i outside of the CMB 
masks. They are used to correct for the incompleteness 
of the survey. The sum is over all pairs whose separa- 
tions are within the angular bins. Correspondingly, the 
weighted cross correlation is given by 
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where, Wj = W(6j) is the weighting function given by 
Eq.@. 

Since we are interested in angular separation of de- 
gree scale (9 = 0(1 )), we can perform the above mea- 
surement for each small patch of the sky, as long as the 
sky patch is much larger than one square degree. 6*df 
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Fig. 3. — The weighted cross correlation of SDSS-DR7 galax- 
ies and WMAP7yr CMB. The direction of the dark flow in the 
weighting is assumed to be I = 282° , b = 22° . The black solid line 
shows result for z < 0.1, and red solid 0.1 < z < 0.2. The black 
dashed line shows the kSZ signal predicted by a single dark flow of 
9530 kms -1 in redshift range z < 0.1, and red dashed line a dark 
flow of 6140 kms - 1 in redshift range 0.1 < z < 0.2. 

is defined as the angle between the center of such sky 
patch and the opposite dark flow direction. For statisti- 
cally isotropic CMB sources (tSZ, the ordinary kSZ, ISW, 
etc.) and cosmic radio and infrared foregrounds, the 
induced cross correlations (contaminations for our pur- 
pose) should be isotropic and be independent of cos #df ■ 
Only the dark flow induced kSZ has the characteristic 
cos 6>df angular dependence. This is the key to separate 
this component from other correlations. 

We estimate the errors of the cross correlation mea- 
surements using the Jackknife method. We divide the 
joint survey region into 1000 subregions which have equal 
effective area (weighted by galaxy completeness). Then 
the Jackknife samples are built by discarding one of the 
subregions each time. The cross correlation is estimated 

by w = wq — (1 — 1/-/V) X)i=i w i: where N is the number 
of subregions, wq is the cross correlation measured from 
the full survey region and lOj from the ith Jackknife sam- 
ple. This estimator is nearly bias free in general cases. 
Especially if w is only a quadratic function as this case, 
the w defined here will be unbiased. The error of this 
estimator is given by 



a 2 (w) 
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where w = YliLi w i/N is the average of iDj. To explore 
the cross talk between different angular bins, we also con- 
struct the full covariance matrix for the Jackknife sam- 
ples, 
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5. MEASUREMENTS AND IMPLICATIONS 



The weighted CMB-galaxy cross correlations are 
shown in figure [3l for two redshift bins z < 0.1 and 
0.1 < z < 0.2. The measurement is pretty noisy. Nev- 
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Fig. 4. — The cos#df dependence of unweighted cross correla- 
tions, w Tg (0|$r>p), f° r redshift range z < 0.1 (upper panel) and 
0.1 < z < 0.2 (lower pannel). The w T 9(e\6 DF ) in the angular 
bins (9) are shown as colored lines. As shown, w Tg (6|0rjp) shows 
a complicated dependence on cosSrjp. However, on sub-regions 
with cos£>df > —0.3 for low redshift bin and cos^df < —0.4 for 
high redshift bin, w T s(6\8 DF ) are nearly linear function of cos $df- 
It is also the case for high redshif sample on sub-regions with 
cos#df > —0.3, except that the linear relation is weak and noisy. 

erthclcss, we detect non-zero correlations at 3.5<7 for low 
redshift bin and 3er for high redshift bin. 

To further quantify the anisotropy of the CMB-galaxy 
cross correlation, we divide the survey region into 10 sub- 
regions according to the value of cos $df and then mea- 
sure w Tg in each sub-regions. The results are shown in 
figure |4j The measured w Tg '(#|#df) shows significant de- 
pendence on #df, which is consistent with a non-zero 
u"w signal. 

As explained earlier, the ensemble average of the 
weighted cross correction vanishes for statistically 
isotropic CMB components and galaxy distributions. 
Hence the detected non-zero cross correlation may imply 
either (1) statistical fluke caused by statistical fluctua- 
tions in CMB and galaxy distribution, (2) anisotropic 
and correlated residual contaminations in CMB and 
galaxy distribution, or (3) the dark flow induced kSZ 
effect that we are searching for. 

5.1. Statistical fluke? 

We have performed a number of tests to check whether 
the detected abnormal anisotropic correlation is caused 
by statistical fluctuations in the intrinsically isotropic 
CMB and galaxy distribution. In particular, we pay at- 
tention to the intrinsic CMB temperature fluctuations, 
the galaxy shot noise and the CMB monopole and dipole. 
These components are uncorrelated with the galaxy dis- 
tribution. However, they do induce statistical fluctua- 
tions to the measured CMB-galaxy cross correlations and 
can in principle mislead the interpretation of the mea- 
surements. 

5.1.1. CMB fluctuations: intrinsic and secondaries 

Firstly we have generated 1000 simulated CMB tem- 
perature maps using the best fit angular power spectrum 
and the beam function provided by the WMAP team. 
Then Gaussian distributed white noise are added to each 
pixel of the map according to the observation times of 



that pixel. We have not included the 1/f instrumental 
noise since it has been removed during map-making pro- 
cedure and it does not correlate with galaxies and hence 
only contributes statistical errors. The weighted cross 
correlation are measured for these simulated maps, and 
the results are shown in the upper panel of figure [5j As 
expected, the intrinsic CMB anisotropies do not bias our 
kSZ signal and the errors are consistent with those of 
Jackknife method. 

The CMB secondaries, such as thermal SZ and 
ISW effect, in principle, are statistically isotropic. So 
under the designed weighting scheme, they contribute 
only the statistical error to our cross correlation mea- 
surements, rather than bias them. Simple order of 
magnitude estimation shows that the induced statistical 
fluctuations are negligible to that of primary CMB. 
The induced statistical error per multiple £ bin is oc 

^J(C? MB + Cf z + C ISW + C> nst + ■ ■ -) w (Cf + Cf' shot )/{2l + 1). 
Here, Ci is the corresponding power spectrum of the 
primary CMB, the thermal SZ effect, the ISW effect, 
instrumental noise, galaxy number density and galaxy 
shot noise, respectively. The subscript "W" denotes the 
corresponding quantity under the anisotropic weighting 
scheme. At angular scales accessible to WMAP, primary 
CMB is the dominant component and hence contributes 
most to the statistical error. So for estimating the sta- 
tistical error of the cross correlation, we can neglect the 
other components. Nevertheless, to fully consider these 
effects one needs to simulate the CMB temperature and 
galaxies simultaneously with cross correlatio ns predicted 
by th e corresponding components. While ICabre et al.l 
(2007) have concluded that within the 10% level, the 
statistical errors coming from simulated CMB maps 
with and without cross correlations with galaxies are 
fully consistent, especially for scales smaller than 10°. 
This result confirms our estimation. So in this work we 
just ignore the CMB secondaries in the statistical error 
estimation. 

5.1.2. Galaxy shot noise 

Secondly, we examine the galaxy shot noise effect by 
generating 100 mock galaxy catalogs using the KIAS- 
VAGC selection function. Each mock catalog has the 
same number of galaxies as the real one. We have mea- 
sured the weighted cross correlation for the mock cata- 
logs and the real CMB data. The results are shown in the 
lower panel of Figure [5j Again, we do not see any bias 
effect and the errors are small enough to be neglected. 
We do not attemp t to include the int rinsic galaxy clus- 
tering effect, since ICabre et al.l (|2007|) have shown that 
using simulated CMB maps one can get sufficient preci- 
sion for this kind of cross correlation measurements with 
low signal-to-noise ratio. 

5.1.3. CMB monopole and dipole 

The CMB dipole are anisotropic over large angular 
scales, which could interfere with our searching for dark 
flow. The CMB monopole has a large amplitude, which 
may induce large statistical fluctuations in the cross cor- 
relation. Hence we remove them in the CMB maps we 
used. Nevertheless, there could exist residual monopole 
and dipole, which may contaminate our cross correlation 
measurement. To test this possibility, we cross correlate 
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Fig. 5. — The systematics coming from intrinsic CMB 
anisotropics (upper) and galaxy shot noise (lower). 
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Fig. 6.— The effect of CMB dipoles and monopoles. The CMB 
dipoles and monopoles are extracted from the real maps using 
the Healpix subroutine "remove_dipole" . The top panel shows the 
weighted cross correlations of galaxies and CMB dipoles, and bot- 
tom CMB dipoles+monopoles. 

the monopole/dipole components of the CMB maps we 
used with the galaxies. The results for the three WMAP 
bands are shown in Figure [H The frequency depen- 
dence presented in the figure comes from the difference of 
dipoles on the maps. For the W band we are focused on, 
the residual dipole contribution to the weighted cross cor- 
relation is less than O.l^tK, so is the residual monopolc. 

5.2. Systematical errors from residual contaminations 

Based on these tests, we conclude that the measured 
anisotropic correlation is unlikely caused by statistical 
fluctuations, such as CMB primaries or the even smaller 
components, TSZ and ISW effects. Rather, it implies 
either a dark flow induced kSZ effect, or caused by the 
residual contaminations in WMAP/SDSS. 

5.2.1. Are CMB intrinsic fluctuations and galaxy 
distribution isotropic? 

One measure of the residual CMB/galaxy contamina- 
tions is the isotropy of CMB and galaxy distribution. 
We measure the sky-dependence of the galaxy and the 
CMB temperature auto-correlation in the same subre- 
gions. The results are shown in figure [JJ None of them 




Fig. 7. — The cos(Sdf) dependence of galaxy and CMB auto- 
correlations. Top: galaxy auto-correlation in redshift range z < 
0.1. Middle: galaxy auto-correlation in redshift reange 0.1 < z < 
0.2. Bottom: CMB auto-correlation. 
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Fig. 8. — The effect of foreground on weighted (upper) and un- 
weight (lower) cross correlation measurement. The black, green 
and red lines show results W, V and Q band respectively. 

show statistically significant sky dependence. This im- 
plies that the residual contaminations are well under con- 
trol and it is unlikely to be the cause of the measured 
anisotropic cross correlation. 

5.2.2. CMB foregrounds 

We compared the cross correlation results using the 
foreground removed CMB maps with the ones which use 
maps without foreground removal, the results are shown 
in figure [3] As can be seen from the figure, the high fre- 
quency band W and V do not suffer much from the fore- 
ground, while the Q band is contaminated badly. After 
the foreground removal procedure, the frequency depen- 
dence is greatly suppressed for both the weighted and un- 
weighted cross correlations, showing that the foreground 
contamination has been largely cleaned, and its contri- 
bution to the cross correlations should be small. 

5.3. The dark flow induced kSZ effect? 

We find that the observed cross correlation can not 
be explained by the reported dark flow scenario either, 
mainly due to two reasons. First, the observed direc- 
tional dependence does not match the prediction of dark 
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flow. As explained earlier, the dark flow induced kSZ 
effect predicts a unique relation w T9 (9\8bf) oc cos6>df- 
However, figure [4] shows that the measured directional 
dependence is rather complicated and can not be well 
described by the above relation. Furthermore, a puzzling 
issue is that over the range cos(#df) > —0.3, the cross 
correlations for both redshift bins do follow the cos(6*df) 
dipole pattern. This is also the case for the 0.1 < z < 0.2 
bin with cos(6*df) < —0.4. Given the noisiness of the 
measurements, we do not attempt to quantify the sta- 
tistical significance of these behaviors. Nevertheless, the 
complicated directional dependence suggests that it un- 
likely, for the cross correlations we found, comes from a 
single dark flow. 

Second, to match the measurements, the required dark 
flow amplitude must be much larger than the ~ 1000 
km/s value reported by IKashlinskv et al.l ([2010). For 
quick comparison, in figure [3] we plot the theoretical 
predictions of kSZ signal coming from dark flow, us- 
ing method described in section [2j We adopt a value 
of 9530 kms -1 over redshift range z < 0.1 (black dashed 
line) and 6140 kms" 1 over redshift range 0.1 < z < 0.2 
(red dashed line) . A dark flow of such large amplitude is 
unrealistic, since it would cause the kSZ signal of clusters 
to be comparable or even dominate over the tSZ signal of 
~ 10 keV clusters and contradict existing SZ observation 
of clusters at low redshifts, especially the measurements 
at the 217 GHz frequency band (e.g. the Planck inter- 
mediate results, Plan ck Collaboration et aLll2012[h 

Hence we rule out the possibility of a dark flow as 
the major source of the observed anomalous anisotropic 
WMAP-SDSS cross correlation. However, the above re- 
sults do not rule out a dark flow of much smaller am- 
pli tude, such as the ~ 10 00 km/s dark flow reported 
by IKashlinskv et al.l (|2010l) . Such dark flow would con- 
tribute a sub-dominant fraction to the measured cross 
correlation and hence do not necessarily contradict with 
the measurements. The low signal-to-noise ratio in our 
measurement prohibits further investigation. The up- 
coming Planck project, with better angular resolution 
and lower instrumental noise, could provide us more in- 
formation on the true origin of the anomalous anisotropic 
correlations. 



6. SUMMARY 

In this paper we report detection of anomalous 
anisotropic CMB-galaxy cross correlations using KIAS- 
VAGC galaxy catalog of SDSS-DR7 and WMAP7 CMB 
data and discuss their implications in the dark flow sce- 
nario. We measure this anisotropy using a specially 
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